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ABSTRACT 


A  Nusselt  type  analysis  was  performed  during  laminar  film  condensa¬ 
tion  on  the  inside  of  a  rotating  truncated  cone.  This  analysis  was 
employed  to  determine  the  condensing  limit  of  a  wickless  heat  pipe, 
rotating  about  its  longitudinal  axis. 

Performance  characteristics  including  the  effects  of  geometry, 
rotational  speed,  and  the  characteristics  of  fluid  are  given.  A 
comparison  is  made  between  the  condensing,  boiling,  sonic,  and  entrain¬ 
ment  limits  for  a  given  heat  pipe  geometry. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  recent  years  considerable  interest  has  been  generated  in  the 
study  of  capillary-pumped  heat  pipes  as  a  means  of  transferring  large 
quantities  of  heat.  At  present,  for  most  pipes  using  ordinary  fluids, 
the  limit  on  the  heat  transfer  appears  to  be  the  rate  that  the  condensate 
can  be  returned  to  the  evaporator  [1],  The  capillary -pumped  heat  pipas 
are  primarily  limited  by  the  small  pressure  differentials  obtainable 
across  the  meniscus  of  the  liquid  in  the  wicking,  which  gives  rise  to 
the  pumping  force.  This  differential  pressure  must  overcome  the  losses 
of  the  liquid  flowing  in  the  wick,  the  losses  in  the  vr.por  flow,  and, 
depending  on  the  orientation,  the  gravitational  forces.  A  second  limit, 
the  onset  of  nucleate  boiling  in  the  evaporator  portion  of  the  wicking 
structure,  may  disrupt  the  capillary  flow  [2),  [3].  Nucleate  boiling  will 
cause  the  formation  of  bubbles  in  the  liquid  condensate.  These  bubbles 
may  then  become  trapped  in  the  wicking  structure,  thus  impeding  the  re¬ 
turn  flow  and  limiting  the  operation  of  the  pipe  (2],  [3]. 

The  limitations  of  the  capillary-pumped  heat  pipes  can  be  overcame 
through  the  use  of  centrlfugsl  acceleration  {&).  A  pipe,  rotating  about 
its  longitudinal  axis,  will  genersta  its  own  gravitational  field,  indepen¬ 
dent  of  orientation,  and  thereby  force  the  condensate  to  return  to  the 
boiler. 

A  rotating  heat  pipe  is  a  closed  hollow  shaft  which  contains  a  work¬ 
ing  fluid.  Rotation  about  its  longitudinal  axis  generates  a  rotational 
force  field.  The  rotational  acceleration,  together  with  a  slight  internal 
wall  taper,  pumps  the  working  fluid  from  the  condenser  to  the  boiler 
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section.  There  is  no  porous  wick  structure  as  In  a  capillary  heat 
pipe.  The  operation  of  a  rotating  heat  pipe  is  shown  schematically 
in  Fig.  1. 

Heat  In  Heat  Out 


Figure  !■  Cutaway  View  of  Typical  Rotating 

Heat  Pipe 

As  the  pipe  is  rotated,  an  annulus  of  working  fluid  is  formed  at 
the  largest  internal  radius  of  the  shaft.  Heat  is  applied  to  this 
portion  of  the  pipe.  As  heat  is  applied,  the  fluid  evaporates  or 
bolls,  forming  a  vapor.  This  vapor  is  then  at  a  higher  pressure  than 
that  at  the  narrower  end  of  the  pipe  (condenser).  This  pressure  differ¬ 
ential  forces  the  vapor  frcn  the  boiler  to  the  condenser  section.  Heat 
is  removed  at  the  condenser  section,  causing  the  vapor  presene  at  that 
section  to  condense  on  the  walls  of  the  shaft.  The  centrifugal  force 
generated  by  the  rotational  motion  of  the  pipe  insures  that  the  con¬ 
densate  will  return  to  the  boiler  eection.  Thus  the  rotating  heat 
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pipe  has  a  closed,  two-phase  heat  transfer  cycle,  and  the  flow  of  the 
condensate  is  not  retarded  by  the  porous  wick  structure. 

The  rotating  heat  pipe  aay  be  applied  wherever  there  is  a  rotating 
part  and  a  need  to  transfer  heat.  A  few  applications  of  the  rotating 
heat  pipe  are  the  cooling  of:  a  high  speed  drill,  the  rollers  for  a 
conveyer  table  of  hot  metals,  the  rotors  of  electric  motors,  and  turbine 
blading. 

B.  THESIS  OBJECTIVES 

The  objectives  of  this  study  were:  a)  to  understand  the  operational 
principles,  b)  to  develop  an  analytical  model  and,  c)  to  predict  the 
performance  characteristics  of  a  rotating  heat  pipe. 
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II.  LIMITS  OF  OPERATION 


As  in  the  capillary  heat  pipes,  the  heat  transfer  limits  of  the 
rotating  heat  pipe  are  controlled  primarily  by  fluid  dynamic  considera¬ 
tions.  [5  j.  The  effect  of  the  following  limits  will  be  discussed:  a) 
the  boiling  limit,  b)  the  entrainment  limit,  c)  the  sonic  limit,  and 
d)  the  condensing  limit. 

A.  BOILING  LIMIT 

When  a  surface  is  heated  in  a  stagnant  pool  of  saturated  liquid, 
three  basic  regimes  of  heat  transfer  are  observed  to  occur.  In  the  low 
range  of  temperature  differentials ,(Tw~T#) , the  heat  flux  is  low  and  the 
principle  mechanism  of  heat  transfer  is  free  convection.  As  the  tempera¬ 
ture  differential  is  increased,  nucleate  boiling  becomes  dominant.  In 
this  region,  bubbles  of  the  vapor  are  formed  on  the  heating  surface, 
grow  in  else,  and  eventually  break  off.  After  breaking  away  they  rise 
to  the  surface  and  burst,  releasing  the  vapor.  The  formation  of  the 
bubbles  and  their  consequent  rising  agitates  the  liquid,  causing  convec¬ 
tive  currents  to  be  set  up.  This  agitation  intensifies  the  heat  transfer 
mechanism  causing  an  upturn  in  the  heat  flux  curve  and  increasing  the 
heat  transfer  coefficient  significantly.  As  the  temperature  differential 
is  Increased  further,  the  bubbles  become  so  numerous,  that  they  begin  to 
coalesce,  and  a  film  of  vapor  can  then  blanket  the  surface.  The  insu¬ 
lating  effect  of  this  vapor  blanket  makes  the  film  boiling  regime  un¬ 
desirable.  A  maximum  heat  flux  exists  for  every  boiler  surface  at  which 
the  transition  occurs  between  nucleate  boiling,  and  film  boiling. 
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During  the  operation  of  the  rotating  heat  pipe,  the  heat  transfer 
coefficient,  as  mentioned  above,  will  be  Increased  due  to  nucleate  boiling. 
A  further  increase  in  the  coefficient  will  be  experienced  due  to  the  ro¬ 
tational  acceleration  of  the  boiler  [6],  [7],  [8].  The  boiling  limit  in 
the  rotating  heat  pipe  will  be  reached  upon  the  departure  from  nucleate 
boiling.  This  will  occur  as  a  sudden  drop  in  the  heat  transfer  coef¬ 
ficient  which  occurs  when  the  maximum  heat  flux  mentioned  above  has  been 
reached.  If  the  power  input  is  held  constant  during  this  time,  the  sur¬ 
face  temperature  of  the  container  will  rapidly  rise,  usually  causing  the 
surface  to  melt.  In  the  literature  this  phenomenon  has  been  termed 
departure  from  nucleate  boiling,  boiling  crisis,  or  burnout  [9].  The 
boiling  limit  of  the  rotating  heat  pipe  may  be  derived  directly  from 
the  work  done  on  conventional  non-rotating  boilers. 

Kutateladze  has  postulated  that  the  transition  from  nucleate  boiling 
to  film  boiling  is  purely  a  hydrodynamic  phenomenon  [10].  Based  on  this 
assumption  he  has  obtained  the  correlation  formula  for  predicting  the 
burnout  heat  flux. 


qt  ,  J/4 

r  =  K  A  g[pf  -pj } 


(i) 


where 


s 

hfB 


=  total  heat  transfer,  BTU/hr 

2 

=  heat  transfer  area  in  the  boiler,  ft 

3 

=  density  of  the  vapor,  lbm/ft 
=  latent  heat  of  vaporisation,  BTU/lbm 

-  surface  tension,  lbf/ft 

2 

-  acceleration  of  gravity,  ft/hr 

3 

-  density  of  the  fluid,  lbm/ft 
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Kutateladze's  experimental  data  Indicates  the  average  value  of  K  is  .14, 
obtaining  from  a  range  of  .13  to  .19  for  various  surface  conditions  of 
horizontal  wires  and  discs  [11].  A  counterflow  configurattion  of  in- 
viscid  vapor  and  liquid  jets  in  the  high  nucleate  boiling  regime  was 
proposed  by  Zuber,  Tribus,  and  Westwater  to  predict  the  burnout  point  [12]. 
Their  results  for  the  burnout  heat  flux  on  a  flat  plate  heater  shows  that 
K  is  equal  to  ^HU  (or. 13).  Experimenters  have  substantiated  Kutateladze's 
formula  and  have  shown  that  the  burnout  heat  flux  Increases  for  rotation 
approximately  as  the  rotational  acceleration  to  the  one  fourth  power  [13], 
[14].  Thus  the  boiling  limit  of  the  rotating  heat  pipe  may  be  adequately 
described  by  the  correlation  formulas  of  Rutateladze  and  Zuber  [9],  [15]. 
The  heat  flux  equation  (1)  may  be  recast  as  a  heat  rate  by  multiplying 
through  by  the  surface  area  of  the  rotating  heat  pipes  boiler.  This  gives 
equation  (2) 

r  r  y* 

V  K^v%AbH  (2) 


where 


K  -  .13  -  .19 


B.  ENTRAINMENT  LIMIT 

During  the  operation  of  the  rotating  heat  pipe,  a  counter  current 
flow  of  the  liquid  condensate  along  the  walla  and  vapor  in  the  central 
core  of  the  pipe  is  set  up.  A  shearing  off  of  the  condensate  by  the 
vapor .holding  up  the  condensate  and  thereby  causing  the  vapor  formation 
in  the  boiler  to  exceed  the  condensate  return,  will  be  defined  as  the 
entrainment  limit.  High  heat  fluxes  during  steady  state  operation  and 
certain  critical  heat  fluxes  during  start  up  conditions  may  possibly 
generate  counter  current  (vapor-condensate)  flow  rates  capable  of 


triggering  entrainment  of  the  condensate.  In  counter  current  flow,  a 
condition  closely  related  to  the  entrainment  problem  has  been  defined  in 
the  literature  as  flooding  [16].  When  this  condition  is  reached,  it  la 
Impossible  to  increase  the  flow  rate  of  the  liquid  component  without 
making  a  corresponding  decrease  in  the  flow  rate  of  the  vapor  component 
and  vice  versa.  Hewitt,  Lacey,  and  Nicholls  have  shown  that  the  flow 
contact  length  ‘is  a  strong  parameter  for  the  flooding  transition  at 
high  liquid  flow  rates.  They  also  have  shown  experimentally  that  as  the 
absolute  pressure  is  increased,  the  flow  rates  attainable  before  the 
flooding  transition  are  higher  [17].  In  order  to  understand  the  state¬ 
ments  and  results  of  Hewitt,  Lacey,  and  Nicholls,  this  discussion  of 
droplet  entrainment  from  a  liquid  film  for  co-current  annular  flow  by 
C.  B.  Wallis  is  presented: 

'Tor  any  given  liquid  flow  rate  the  gas-liquid  Interface  is 
characterised  by  the  following  physical  appearance.  At  aero  gaa 
velocity  the  liquid  has  a  smooth,  mirror-like  appearance  as  it 
moves  along  the  duct.  At  a  low  gas  velocity  the  surface  etarts 
to  break  up  into  a  ripple  pattern.  As  the  gas  velocity  is 
further  increased,  the  ripples  start  to  bunch  up  and  a  wave  is 
produced  which  appears  to  move  at  a  velocity  several  times  that 
of  the  entrained  film.  This  wave,  referred  to  as  a  roll  wave, 
moves  down  the  duct  on  top  of  the  entrained  film.  At  still 
higher  gas  velocities  these  roll  waves  grow  smaller  (yet  au>re 
frequent)  and  droplet  entrainment  becomes  evident.  These  waves 
are  a  precondition  for  the  onset  of  entralnsMnt  and  the  source 
of  the  droplets.  The  wave  tops  are  sheared  off  by  the  gas 
flowing  over  them.  The  roll  waves  appear  to  *roll'  in  the 
direction  of  the  gas  flow."  [18] 


The  rotational  acceleration  in  the  rotating  heat  pipe  will  be  a 
significant  parameter  in  suppressing  the  formation  of  a  rolling  wave  aa 
described  by  Wallis  [18].  It  should  act  in  the  same  manner  as  the  in¬ 
creased  absolute  pressure  delayed  the  onset  of  flooding  in  the  experiment a 
of  Hewitt,  Lacey,  end  Nicholls  1)7]  In  addition,  since  the  rotating  heat 
pipe  is  a  sealed,  self  contained  unit,  the  internal  pressure  may  be 
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regulated  to  such  a  value  as  to  ensure  no  liquid  entrainment  will 
occur. 

At  the  present  time  there  is  no  satisfactory  theory  available  to 
predict  the  onset  of  entrainment,  but  the  experimental  data  of  Hewitt, 
et  al,  can  be  used  to  predict  the  entrainment  limit  [17]. 


C.  SONIC  LIMIT 

As  the  heat  flux  in  the  rotating  heat  pipe  is  increased,  it  is 
possible  that  the  flow  rate  of  the  vapor  down  the  container  may  be 
limited  by  the  attainment  of  a  choked  flow  condition  for  the  vapor. 

This  choking  of  the  vapor  flow  will  definitely  limit  the  amount  of 
energy  the  vapor  can  transport  down  the  pipe,  thus  limiting  the  opera¬ 
tion  of  the  pipe.  The  heat  transfer  mechanism  employed  in  the  rotating 
heat  pipe  depends  on  the  latent  energy  of  the  working  fluid.  The  heat 
transfer  rate  in  the  pipe  may  then  be  expressed  as 

ql  =  "v  hfg  O) 


where 


rhy  «  mass  flow  rate  of  the  vapor  (Ibm/hr) 

Pram  continuity,  the  mass  flow  rate  of  the  vapor  may  be  expressed  by 


«\=fvuvA 


(4) 


where 

uy  •  velocity  of  the  vapor  (ft/sec),  and 

2 

A  •  cross  sectional  area  for  flow  (ft  ) 

Combining  equations  (3)  and  (4)  yields 

q.= O) 
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If  Che  cross  sectional  area  of  Che  pipe  is  known  where  Che  choked  condi¬ 
tion  occurs,*  the  sonic  limi cation  on  Che  heat  rate  nay  Chen  be  com- 
puted  using  equation  (5)  with  the  vapor  velocity  taken  as  Che  sonic 
velocity. 

D.  CONDENSING  LIMIT 

The  condenser  nay  limit  the  rate  of  return  of  the  condensate  to  t&e 
boiler  by  two  mechanism;  a)  the  rate  that  energy  My  be  extracted  from 
the  system,  and  b)  the  maximum  flow  rate  of  the  condensate.  At  the 
present  no  solution  exists  in  the  literature  for  condensation  on  a 
truncated  cone.  The  solution  for  condensation  on  a  rotating  cone  has 
been  carried  out  by  Sparrow  and  Hartnett  [19].  Their  solution  employs 
a  boundary  layer  similarity  approach.  As  is  stated  by  Sparrow  and 
Hartnett,  the  solution  applies  only  to  cones  that  are  not  too  slender 
and  in  a  region  sufficiently  removed  from  the  apex  of  the  cone.  Ae  can 
be  seen  in  Fig.  1,  the  rotating  heat  pipe  will  generally  employ  a 
slender,  truncated  cone  for  a  condenser  section.  If  Sparrow  and  Hartnett's 
solution  for  condensation  on  a  cone  were  employed,  then  the  boundary  condi¬ 
tion  of  the  rotating  heat  pipe  would  not  be  fulfilled,  namely  that  the 
velocity  of  the  condensete  is  sero  at  the  narrowest  section  of  the  con¬ 
denser.  Hence,  a  truncated  cone  solution  Is  necessary.  However,  since 
a  similarity  solution  for  the  truncated  cone  would  not  apply  in  the 
region  near  the  narrow  end  of  the  condenser,  due  to  the  boundary  layer 
assumptions,  the  following  Nusaelt-type  analysis  was  Mda. 


As  slender  pipes  are  anticipated  (smII  cone  angles)  then  as  a 
good  first  approximation,  it  can  be  assumed  that  choked  flow  occurs  at 
the  exit  area  of  the  boiler. 
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III.  THEORETICAL  ANALYSIS  OF  FILM  CONDENSATION 
OH  THE  INSIDE  OF  A  ROTATING -TRUNCATED  CONE 


In  attacking  the  problea  of  condensation  on  the  inside  of  a  trun¬ 
cated  cone  the  following  assumptions  are  made. 

A.  ASSUMPTIONS 

a)  Film  condensation,  not  dropwise  condensation,  occurs  in  the 
condenser  section.  Thus  the  condensing  liquid  wets  the  container  sur¬ 
face  and  condensation  occurs  at  a  sufficient  rate  to  allow  a  continuous 
film  over  the  surface.  The  assumption  of  filmwise  condensation  is  felt 
to  be  the  conservative  choice  as  dropwise  condensation  has  been  shown 
to  give  higher  heat  transfer  rates  (20). 

b)  The  film  of  condensate  undergoes  laminar  r 1 ow .  Cannon  and 
Kays,  using  a  pipe  rotating  about  its  longitudinal  axis  with  annular 
flow,  have  reported, 

"Tube  rotation  apparently  has  the  effect  of  stabilising  a 
laminar  flow  so  that  a  higher  through -flow  Reynolds  number  is 
obtained  before  transition  occurs."  (21]. 

c)  The  fluid  properties  are  constant.  Heat  pipes  are  nonsally 
operated  as  nearly  isothermal  devices.  This  condition  should  be  easily 
achieved  in  the  rotating  system  and  thus  the  fluid  properties  will  very 
nearly  be  a  constant. 

d)  The  subcool ina  of  the  condensate  may  be  neglected .  A  small 
correction  factor  is  normally  applied  to  the  latent  heat  of  vaporisation 
term  In  order  to  take  into  account  the  effect  of  subcooling.  The  neces¬ 
sity  for  this  correction  is  determined  by  the  comparison  of  the  term 

C  AT  (C  Is  the  specific  heat  at  constant  pressure,BTU/lbm  deg  F,  and 
P  P 

AT  is  the  difference  between  the  wall  temperature  and  the  saturation 
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temperature,  deg  F)  with  the  latent  heat  of  vaporitation.  Normally 

C  AT  le  orach  leae  than  h.  and  thua  may  be  neglected, 
p  fg 

e)  The  momentum  change a  through  the  condenaate  are  email  and 
thua  there  is  eaaentlally  a  etatie  balance  of  forcea.  Sparrow  and  Gragg 
have  shown  that  the  momentum  effects  are  negligible  for  ordinary  fluids 
in  laminar  film  condensation  on  a  flat  plate  [22).  Assumption  1),  a  thin 
film,  also  indicates  that  momentum  changes  may  be  neglected. 

f)  The  vapor  exerts  no  drag  on  the  motion  of  the  condensate; 
there  is  no  interfacial  shear.  Experiments  have  shown  that  for  Prandtl 
numbers  of  ten  and  higher  the  effect  on  heat  transfer  of  interfacial 
shear  is  negligible.  For  Prandtl  numbers  near  one  the  effect  la  only 
slight  (within  five  percent)  [23]. 

g)  Thera  is  a  linear  distribution  of  temperature  through  the 
condensate  film.  Rohaenow  has  shown  that  for  Prandtl  numbers  above  one, 
the  non-linear  distribution  of  temperature  through  a  thin  film,  due  to 
energy  convection,  is  hardly  distinguishable  [24].  Thus  energy  convec¬ 
tion  la  oeglected  and  heat  transfer  occurs  by  conduction  only. 

h)  The  centrifugal  force  ia  much  greater  chan  the  force  of 
gravity  and  thua  gravity  nay  be  neglected. 

1)  The  thickness  of  the  film  is  much  less  than  the  radius  of 
curvature  of  the  condenser  wall. 

))  The  vapor  apace  is  essentially  at  one  pressure.  P^. 

There  la  no  radial  distribution  of  pressure  in  the  vapor  space  and  the 
axial  distribution  of  pressure  is  nearly  a  constant  [25],  [26]. 

k)  The  density  of  the  fluid  is  much  greater  than  the  density  of 
the  vapor. 

l)  The  rotating  heat  pipe  is  operating  at  a  steady  state  condi¬ 


tion. 


■)  Velocity  gradients  in  the  circuaferential  direction  rela¬ 
tive  to  the  pipe  vail  are  negligible. 

By  using  the  above  assumptions  and  the  coordinate  system  shown  in 
Fig.  2,  an  analysis  similar  to  Nusselt's  original  film  condensation 
theory  may  be  used  [27]. 


Figure  2  Coord'note  $y*t«m  0nd  Geometry 
for  Condtnting  Limit 


A  force  balance  may  be  taken  on  an  infinitesimal  fluid  element 
in  the  X -direction  as  shown  in  Fig.  3.  this  leads  to  an  expression 
for  the  momentum  equation  within  the  condensate. 
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Figjre  3  Force  Boloncein  X- Dire  cl  ion  on  on 
Infimtesimol  Fluid  Element 


1.  HOHBWnJH  BWAT10M  (X-dir«ction) 

lF,«2»rP<ly-2»r  {p,41a,j,(-r2.rdi  + 

{r  *  i/}Zwr  i’  *  fifw1  2  wrl  it  ty  tin  f 

.where 

t  •  pressure,  (lbf/ft2) 
r  •  radius, (ft) 

X  •  co-ordinate  Maturing  distance  along  the  surface ,  (ft) 

T  •  co-ordinate  Maturing  distance  normal  to  surface,  (ft) 
T  "  Shear  stress,  (lbf/ft2) 


(6) 


n 


<fc  ”  half  cone  angle  (rad) 

,  2 

pf  *  density  of  fluid,  (Ibo/ft  ),  and 
w  ■  angular  velocity,  (1/hr) 


This  equation  reduces  to 


0s 


dr 

dy 


dP 


+  /3ju7f  Slf'^ 


C.  MOKEHTTUM  SQUATa-ON  (T -direction) 

In  like  manner,  using  Fig,  4,  a  force  balance  in  the 
neglecting  the  centrifugal  head  due  to  the  vapor  since 

£F  =27rrPdx  -  2rf  p+  ~  dyj dx -p^xZir  cos f  dxdy 


This  rsducee  to 


(7) 

Y -direct! on, 
>:>PV  »  give*: 

(8) 

(9) 


Making  use  of  Fig.  2,  the  radius  at  any  point  in  the  fluid 


condensate  may  be  expressed  as 


r  (;  y)  =  Rp  +  X  sin  ^5  -  Y  cos  <p 


(10) 


where 


R  *  minimum  wall  radius  in  the  condenser  section  (ft) 
0 


Substituting  equation  (10)  into  (9)  for  r(x,y)  and  integrating  gives 

&  p  8 

fjr  dy  -- pLucosj)  f(R  +  Xsin<£  -  Y  cos<p  dy  (ll) 

P  y  Y 


where 


*  pressure  of  the  vapor  (lbf/ft  ) 


8  *  film  thickness  (ft) 


This  yields 


p=  {[r0+ Xsin^>][s-Y]-  ~  cn*rf<|V-  Y2]}  (12) 


Now  differentiating  equation  (12)  with  respect  to  X  gives 


-37-  s  +AuT cc  I  S0~- —  +  8smA+Xb»n^~^  (13) 

dx  dx  rf  ^  1  °d  x  T  dx 

-  Y  sin<$>  -  8  cos  £  ^7*} 

»p  r  d  x  j 

If  the  result  of  equation  (13),  * —  ,  is  substituted  into  equation 

0  x 

(7)  it  becomes 


dPu 


0=  ~  "dxA,~/5fa'2cos^{‘^dx  +^Sin^  +  Xsin  P  Z~  '  Ysm^-Scos^J^-}  (14) 


f)  -‘^.r  R0+  as-'  ‘  -  Ycos  /•)  +-~ 
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Integrating  equation  (14)  over  the  limits  from  Y  to  8  and  r  to  0, 
then  rearranging  yields 


I  .  ’  *  -s 


-8cosf  +  PjU  sin  / iR0+ Xsi .•  -y} 


where  in  equation  (IS)  for  laminar  flow 

r  *'  t‘i  r; 


u^  ■  velocity  of  the  fluid,  (ft/hr) 


fx f  m  viscosity  of  the  fluid,  (lbm/ft  hr) 

Equation  (15)  may  now  be  solved  for  the  velocity  profile,  ,  by 
integration  over  the  limits  0  to  Y  and  0  to  , 

u  -  _  J_s1£v%y-  S^{ +Ssin^+xsin<^~  (] 

f  Pfdx  1  2  )  Pj  1  bci*  'f  r  dx  v 

£  ){*Y- f }  ♦  ^{f0+  XS,n^}{SY-^2} 


Applying  assumption  j),  (  vapor  space  at  one  pressure,  ),  allows 

equation  (16)  to  be  rewritten  as 


Uf  =  {^Y -  £  }  { F  i;  /- -  cos  /  ~]  +  X sin -/>!* -in  / -cos^  4^-]  (17) 

-  8c  )S1  jsin  £  -cos^^-J} 


Note  that  in  equation  (17)  u  »  u^x.y)  since?  »  8  (x). 
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D.  CONTINUITY  EQUATION 


E.  ENERGY  EQUATION 

If  the  film  surface  temperature  is  T  and  the  wall  is  at  T  ,  then 

0  W 

the  heat  transfer  by  conduction  to  an  element  of  surface  area  dA&  is 


dq,  = 


{R0t  X sin  f  }  In  { R)>» xsto$-! f  cos<f ) 


(22) 


where 


dA£  «  ?.  it }  Kq  <  X  sir  /■}  lx 

A  ■  heat  transfer  area  in  the  condenser,  (ft) 
c 

kf  *  thermal  conductivity  of  the  condensate  film, 
(BTU/hr  ft  deg  P) 

T  *  saturation  temperature,  (deg  F) 

8 

*  inside  wall  temperature,  (deg  F) 


Therefore 


‘tqt* 


which  can  be  written  as 


2vKfCOSf  lTs~V 


U 


'"fid 


%±1±l4 _ _) 

Xsin  /  -  ^  r>s  f  ' 


(23) 


ln  i1 " 

The  assumption  of  a  thin  film,  assumption  i),  gives 

£  r  /• 

K  <  :• 


(24) 


Therefore 


In  1 1  — 


ft I  ) 
R  +  X  si  n  /  ' 


,S  r  <;  <f\ 
*■  X  sm  1 
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final  raaulta  ara 


2 

The  local  heat  transfer  coefficient,  h,  (BTU/hr  ft  deg  F) ,  it  given 

by 


G.  MEAN  HEAT  TRANSFER  COEFFICIENT 

The  Man  heat  transfer  coefficient  over  the  whole  surface  area 
of  the  truncated  cone  is  given  by 


h  = 
m 


t-C 

/  2w,\ 


,t  / ,  • i . 


(32) 


where 

1^.  -  length  along  the  wall  of  the  condenaer,  (ft) 
Integrating  equation  (32)  (see  Appendix  B)  yields 


m 


I  f  2  r.  . _ iM 

{2fV  +  Lrsin  +}  i  ? 


(33) 


* '  ^ .  4  i.  i  •  >•' 
r 


V4 
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IV.  DETERMINATION  OF  THE  CONDENSING  LIMIT 


The  analysis  of  fils  condensation  on  a  truncated  cone  say  be  ap¬ 
plied  to  the  condenser  section  of  the  rotating  heat  pipe.  This  then 
will  give  an  equation  which  will  predict  the  condensing  lisit  for  the 
rotating  heat  pipe.  The  total  heat  rate  in  the  condenser,  q£,  is  given 
by 

=  hmir  4  Y  V  i2  V  (34) 


The  value  of  h  is  now  substituted  into  equation  (34)  and 
0 


R0+Lcsm 


>3  6T3  » 

*]  -Ro* 


This  then  is  the  operating  equation  for  predicting  the  Halting 
heat  rate  for  the  condenser.  In  the  equation,  however,  the  quantity 
Tv  is  not  a  readily  controllable  quantity.  The  inside  wall  taaperature 
will  actually  adjust  to  a  taaperature  corresponding  to  the  heat  rate 
(the  saturation  taaperature)  and  the  outside  wall  taaperature.  As  the 
heat  rate  in  the  rotating  h^at  pipe  is  increased,  the  Internal  pres¬ 
sure  and  saturation  temperature  will  increase,  accordingly  requiring 
that  the  inside  wall  temperature  adjust  for  this  Increase.  Thus,  the 
outside  wall  temperature  of  the  condenser  will  be  the  controllable  and 
measurable  quantity.  The  inside  wall  temperature,  T^,  then  suiy  be  found 
as  a  function  of  the  outside  wall  temperature,  T  ,  and  the  heat  rate 

VO 

by  use  of  Fourier's  conduction  law  through  the  consenser  wall, 

dV  {"2"{vx',n*}  {VT„<}d* 


(34) 


where 


■  thermal  conductivity  of  the  well,  (BTU/hr  ft  deg  F) 

^t  •  condenser  well  thickness,  ft 

Once  egein  the  cylindrical  effects  have  been  neglected  due  to  the  thin 
well  of  thickness  t.  Integrating  equation  (36)  from  0  to  qfc  and  0  to  L£ 
yields 

V  ^2irLttVU{R'  +  fS,n*>  (37) 


A.  OPERATING  EQUATIONS 

The  operating  equation  for  the  condensing  limit  is  then  an  itera¬ 
tive  solution  between  equation  (35)  and  (37). 


3  I 

*,P 


iw]' 


sir f  / 


1/4  8/3  3/3  . 


(35) 


*"*}  (37) 


At  this  point  it  should  be  noted  that  although  operating  equa¬ 
tions  for  the  condensing  limit  of  the  pipe  Itself  have  been  derived, 
as  an  engineering  consideration,  the  outside  cooling  mechanism  for 
the  pipe  is  extremely  important.  Indeed,  it  may  became  the  limiting 
consideration.  This  then  gives  a  third  equation  to  be  satisfied.* 

Of :  V  "4  iT«o~T‘h  (»> 


*This  is  seeming  that  the  tx  *  and  well  thicknesses  are  very  thin, 
and  thus  the  areas  are  approximately  equal. 
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where 


2 

h  -  coolant  heat  tranafer  coefficient  (BTD/hr  ft  deg  F) 
o 

Tf  *  coolant  teaperature,  (deg  F) 


V.  PREDICTED  PERFORMANCE  CHARACTERISTICS 


In  the  actual  application  of  tha  rotating  haat  pi pa,  varloua  aathoda 
of  cooling  tha  condanaar  faction  arc  anticipated.  Tha  aaaunption  la 
awde  her* , however,  that  tha  cooling  »a chan lan  employed  la  capable  of 
aattlng  tha  outalda  aurfaca  tanparatura  of  tha  condanaar  at  a  con- 

atant  value,  T  Equation  (38)  la  thua  a 1 Inina tad  fren  tNa  aolutlon. 

vo 

Rearranging  tha  regaining  condanalng  Unit  aquatlona  givaa 

%  3 


(35) 


V  2ir 


(37) 


Thaaa  aquatlona  are  a  function  of  tha  phyalcal  propartlaa  of  tha  working 
fluid,  tha  thermal  conductivity  and  the  geometry  of  tha  pipe,  and  the 
operating  condition*  of  rotational  angular  velocity  and  outalda  vail 
tanparatura  T  ■ 

W 

A  rotating  haat  pipe  with  the  phyalcal  diaenaiona  aa  ahovn  In  Pig. 

3  vaa  choeen  In  order  to  danonatrate  tha  predicted  performance  of  a 
given  pipe.  Tha  boiler  eectlon  vaa  aodeled  aa  ahovn  in  Pig.  5  in  order 
to  inaura  that  an  adequate  aupply  of  working  fluid  would  be  available 
for  an  actual  pipe.  Alao  with  thla  nodal  the  rotating  haat  pipe  boiler 
la  cloaaly  approximated  by  tha  rotating  boiler  taatad  at  Levi*  teaaarch 
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Center  (7].  The  reeulte  of  the  Lewie  teste  can  be  used  to  design 
and  cheek  the  operation  of  the  rotating  heat  pipe  boiler. 


Figure  5'  simp  lified  Model  of  o  Rototmg 


Heot  Pipe 


The  outside  wall  temperature  of  the  condenser,  T  was  chosen  as 

WO 

70  degrees  Fahrenheit.  For  the  saturation  temperature  of  a  parti¬ 
cular  fluid,  the  fluid  properties  were  determined.  An  IBM  Dual  360 
digital  computer  was  employed  In  Iteratively  solving  the  condensing 
limit  equations  (33)  and  (37).  For  a  given  angle,  f,  and  rotational 
speed,  ,  the  Inside  wall  temperature,  T^,  was  varied  until  the 
heat  rates  of  equations  (35)  and  (37)  sgreed  within  10  ITU  per  hour. 
It  should  be  noted  that  the  condenser  length,  10  Inches,  and  the 
radius  of  the  boiler  exit,  1.25  Inches,  were  held  constant  while  the 
eagle  #  was  varied. 
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A.  EFFECT  OF  FLUIDS 

The  choice  of  the  working  fluid  in  a  rotating  heat  pipe  la  dependent 
upon  several  physical  properties  of  the  fluid  as  veil  as  the  operating 
conditions  desired.  Various  working  fluids  may  be  compared  by  the  use 
of  a  dimensional  parameter,  N,  defined  as  "the  figure  of  merit." 

,  .3  p2.  JIZ 

N  =  {  ^  fg}  (39) 

The  heat  rate  will  be  proportional  to  the  figure  of  merit  of 
various  fluids  for  a  given  saturation  temperature.  Table  I  gives  the 
variation  in  N  for  three  fluids  at  a  saturation  temperature  of  122 
degrees  Fahrenheit. 

TABLE  I 

Figure  of  Merit 


Fluid  N 

Water  20 

Ethyl  Alcohol  5 

Freon  12  2.4 


It  should  be  noted  that  for  ordinary  fluids,  water  is  the  most 
desirable  working  fluid.  This  is  primarily  due  to  the  combination  of 
a  relatively  high  thermal  conductivity  (which  is  cubed)  coupled  with  a 
high  latent  heat  of  vaporisation. 

The  vaporisation  curve  of  a  given  fluid  will  determine  the  operating 
temperature  range  of  the  fluid.  The  heat  rate  attainable  will  then  be 
determined  by  the  choice  of  fluid  and  the  operating  temperature  of  the 
pipe  as  shown  in  Fig.  6.  It  should  be  noted  that  in  Fig.  6  the  geometry 
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of  the  rotating  heet  pipe  1*  as  shown  in  Fig.  5  with  a  constant  rotation¬ 
al  apead  of  2700  RPM,  and  the  half  cone  angle ,  0,  of  3  degrees. 

B.  EFFECT  OF  GEOMETRY 

For  the  given  geometric  configuration  shown  in  Fig.  5  the  half  cone 
angle,  0,  was  allowed  to  vary  for  a  given  operating  condition  and  work¬ 
ing  fluid.  Figure  7  shows  a  plot  of  the  heat  rate  versus  the  half  cone 
angle  for  water  at  a  saturation  temperature  of  212  degrees  Fahrenheit. 

The  rotational  speed  was  allowed  to  vary  from  300  to  3000  RPM. 

Figure  7  shows  that  a  maximum  heat  rate  is  achieved  at  a  particu¬ 
lar  half  cone  angle  for  a  given  rotational  speed  (the  dotted  line). 

This  peaking  is  due  to  the  model  used.  The  boiler  exit  diameter  was  held 

constant  while  the  half  cone  angle,  0,  was  varied.  Since  L  .cos  0,  the 

c 

condenser  axial  length  was  also  held  constant,  the  minimum  radius  of  the 

condenser,  R  .  varied.  As  0  increased  R  decreased  and  thus  the  heat 
o  o 

rate  versus  the  half  cone  angle  curve  has  a  maximum  point. 

Figure  7  illustrates  that  for  the  model  chosen,  as  the  rotational 
speed  of  the  pipe  is  Increased,  the  half  cone  angle  necessary  to  main¬ 
tain  the  maximum  heat  rate  decreases. 

C.  EFFECT  OF  ROTATIONAL  SPEED 

The  rotational  speed  versus  the  heat  rate  is  shown  in  Fig.  8  for 
various  half  cone  angles,  0.  This  curve  is  presented  for  water  at  212 
degrees  Fahrenheit  in  order  to  determine  the  most  desirable  operating 
speed  range.  One  of  the  assumptions  for  the  condensation  on  a  truncated 
cone  was  that  the  centrifugal  force  is  much  greater  than  the  gravity 
force  (assumption  h).  For  the  model  shown  in  Fig.  5,  this  condition  is 
estimated  to  occur  at  approximately  700  RPM. 
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Heat  Transfer  Rote 


...  _  Condensing  Limit 

rigure  ('  Heat  Transfer  Rate  vs  Half  Cone  Angle 
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Figure  8-  Heat  Transfer  Rote  vs  Rotational 
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D.  OVERALL  CAPABILITY 

A  rotating  heat  pipe  with  the  geometry  shown  *n  Pig.  5,  operating 
at  2700  RFM,  with  a  half  cone  angle  of  1  degree,  and  using  water  as  the 
working  fluid  was  chosen  to  demonstrate  the  effects  of  the  boiling,  en¬ 
trainment,  sonic,  and  condensing  limits  on  the  maximum  heat  rates  at¬ 
tainable. 

The  boiling  limit  curve  was  determined  from  equation  (2),  using 
both  Kutateladze's  value  of  K  equal  to  .19  and  Zuber's  value  of  w  /24 
(.13),  thus  giving  a  range  of  heat  rates  for  the  boiling  limit.  Since 
the  boiler  surface  of  the  rotating  heat  pipe  may  be  approximated  by  a 
flat  plate,  for  large  radii  of  curvature,  Zuber's  predicted  curve  should 
be  the  most  applicable  for  the  model  chosen. 

As  mentioned  earlier,  no  suitable  theory  is  available  to  predict  the 
entrainment  limit,  so  that  the  experimental  results  of  Hewitt,  Lacey, 
and  Nicholls  have  been  used  { 17  ] .  Their  data  was  taken  from  a  1.25  inch 
diameter  stationary  tube  in  annular  flow  with  a  counter-current  air- 
water  contact  length  of  9  Inches.  They  present  flooding  data  for  the 
vapor  mass  flow  rate  as  a  function  of  the  liquid  mass  flow  rate.  The 
assumption  of  steady  state  operation  in  the  rotating  heat  pipe  allows 
the  Hewitt,  Lacey,  and  Nicholls  data  to  be  entered  at  the  condition  of 
equal  liquid  and  vapor  mass  flow  rates.  The  mass  flow  rate  thus  obtain¬ 
ed  was  then  adjusted  to  agree  with  an  average  rotating  heat  pips  dia¬ 
meter  of  2  inches  and  a  counter -current  steam-water  flew  condition. 

This  adjusted  mass  flow  rate  of  vapor  was  then  used  in  equation  (3) 
to  give  the  entrainment  limit. 

(3) 
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Figure  9  Heat  Tranefer  Rate  v*  Saturation  Temperature 


Equation  (5)  was  used  in  computing  the  sonic  limit.  The  asswption 
was  made  that  choking  occurred  at  the  exit  of  the  boiler  section  for  a 
given  saturation  temperature  in  the  condenser. 

The  condensing  limit  was  obtained  by  a  solution  of  equations  (35) 
and  (37)  as  described  on  page  33. 

Figure  9  gives  the  plot  of  the  rotational  heat  pipes  heat  rate  as 
a  function  of  the  saturation  temperature  for  die  limits  Just  discussed. 

For  each  of  the  limit  curves  shown  in  Fig.  9  there  is  an  associated 
uncertainty.  Even  with  this  consideration,  it  is  evident  that  the  ro¬ 
tating  heat  pipe  shown  in  Fig.  5  will  be  condenser  limited,  for  a  given 
saturation  temperature.  The  upper  limit  on  the  heat  rate  attainable 
for  .»  certain  working  fluid  will  generally*  be  governed  by  the  ability 
of  the  container  to  withstand  the  saturation  pressure. 


It  is  possible  that  a  container  could  withstand  a  certain  satura¬ 
tion  pressure,  and  yet  not  be  compatible  with  the  associated  wall 
temperature. 
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VI.  CONCLUSIONS 


1.  The  fila  thickness  end  local  heat  transfer  coefficients  for 
condensation  on  a  rotating,  truncated  cone  aay  be  determined  by  a 
Nusselt  type  of  analysis.  The  resulting  fila  thickness  and  local  heat 
transfer  coefficient  are  not  constant  as  for  the  rotating  cone  and  disc, 
but  vary  with  axial  position. 

2.  The  predicted  perfonsance  for  the  heat  pipe  geometry  chosen  in 
this  study,  and  operating  with  water  as  the  working  fluid,  is  con¬ 
denser  Halted,  for  a  different  geometry,  the  entrainment  Halt  aay 
became  the  dominant  Halt.  The  sonic  and  boiling  Halts  should  pre¬ 
sent  no  practical  limitation  on  the  rotating  heat  pipe. 

3.  Of  the  ordinary  fluids,  water  has  the  best  characteristics 
for  use  as  the  working  fluid. 

4.  The  actual  operation  of  the  rotating  heat  pipe  aay  be  Halted 
by  the  thermal  conductance  of  the  outside  surface,  (l.e.,  upon  the 
method  of  cooling  the  rotating  heat  pipe). 
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VII.  RECOfCHMTIOWS  FOR  FP1THKR  STOPT 


A.  ANALYTICAL 

1.  An  optimization  gtudy  should  b«  conductsd  to  datamlita  tha 
effect  of  the  internal  geometry  (boiler  exit  area  and  condenser  shape 
as  well  as  length)  on  the  rotating  heat  pipe  condensing  llnlt. 

2.  The  restrictions  of  negligible  Interfacial  shear,  negligible 
■omentum  changes  in  the  condensate,  and  constant  vapor  pressure  should 
be  removed  from  the  present  analysis.  This  then  would  require  a  solu¬ 
tion  of  the  continuity,  momentum,  and  energy  equations  for  the  con¬ 
densate  as  well  as  the  vapor.  The  equations  of  the  condensate  and  the 
vapor  would  be  coupled  through  the  boundary  condition  that  the  shear  at 
the  Interface,  in  each  phase,  must  be  equal. 

The  swirl  component  of  the  flow  of  condensate  should  also  be  taken 
into  account.  This  would  then  require  that  the  problem  be  solved  in 
three-dimensions,  probably  through  the  use  of  spherical  co-ordinates. 

3.  The  exact  location  of  the  choked  flow  condition  in  the  rota¬ 
ting  heat  pipe  should  be  determined  in  order  to  more  accurately  chart 
the  sonic  limit. 

4.  A  study  should  be  undertaken  to  determine  accurately  the  en¬ 
trainment  limit  of  the  rotating  heat  pipe  for  the  condenser  geometry 
shown  in  Pig.  3. 

B.  DCPEtDCDCTAL 

1.  An  experimental  model  should  be  built  and  the  condensing  heat 
rate  limit  tv-rre  for  the  pipe  should  be  experimentally  verified  for 
various  working  fluids. 
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2.  A  comparison  of  filawise  and  dropwise  condensation  in  the 
condenser  section  should  be  made. 

3.  A  further  extension  of  the  experimental  procedure  should  be 
to  have  a  condenser  wall  with  axial  grooves  in  it.  This  then  would 
test  the  contributions  of  interfacial  shear  and  swirl  on  the  condensing 
limit. 
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APPENDIX  A 


SOLUTION  FOR  THE  FILM  THICKNESS . 8 (X) 


Equation  (30)  nay  be  written  in  the  foza 


T,J 

p,  tv  hf  sin  ^ 


3  {rq*  X  sin<^|  5  +2  8  s»n<£ 


(A-I) 


Let 

K  = 

R  =  R  +\  sm  >f. 

U  =  84 

Therefore 


4A^in* 


±  K 
3  sin 


}  dR  =  0 


Now  let 


Thle  gives 


K 

sin  $ 


_3  dS  dR 

9  ~b~  R 


Integrating  this  gives 


(A-2) 


(AO) 


(  A-4) 
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where  C  is  an  arbitrary  constant.  Transforming  equation  (A-4)  back, 
into  the  original  dimensions 


f  e_  £  ±  k  r  c _ 

l  3  3  "sm$J  ‘  R0  +  x  s.n  <5b 


(A -5) 


In  owder  to  determine  the  arbitrary  constant,  C,  the  boundary  condi¬ 
tion  8  -  0  at  X  -  0  will  be  used.  The  constant,  C,  then  is 


± 

3 


K 

sin  <£ 


3/e 

v 


(A-6) 


Now  substituting  the  original  value  of  X  and  the  value  of  C  into  equa¬ 
tion  (A-5),  then  rearranging,  the  final  solution  for  8  Is 


8  = 


4  —  \ ;  4--^ — 

^  p*  u/sin  <£>  i  -  RQ+Xsm  <f> 


(A-7) 


The  solution  for  condensation  on  a  rotating  cone,  Sparrow  and 

Hartnett's  solution,  can  be  easily  obtained  by  the  use  of  the  boundary 

condition  that  8  ■  0  at  X  *  0  for  R  ■  0.  This  condition  would  give, 

o 

from  equation  (A-5),  that  C  *  0.  The  solution  for  delta  is  then  the 
solution  obtained  by  Sparrow  and  Hartnett. 
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APPENDIX  B 


SOLUTION  FOR  THE  MEAN  HEAT 
TRANSFER  COEFFICIENT*  h™ 


Equation  (33)  may  be  written  In  the  fora 


Working  with  juat  the  integral  portion  of  equation  (B-l) 


_  V,  k-xwf  r 

z  1  '  ~T'rT 

i  Rq*  X  s  in  <p 

let 


U«  Rc  +  Xsm  / 


o.-> 


Equation  ( B-2)  then  becomes 


let 

Therefore 


It  should  be  noted  that 


(B-3) 


<»-4) 
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Carrying  out  the  integration  therefore  yields 


,  I  r  9  f 

Z=2^tV,A} 


(B-5) 


Finally 


uSinf 


{[Rc+Lcsm^  - 


(B-6) 


Substituting  2  back  into  equation  (B-l),  the  average  heat  transfer  co¬ 
efficient,  h  ,  is  given  by 

ID 


2  ktptui  hx 


1 _ r  £  n<g  i 

{2^+LjSin^}  13  Mf  [TS“TJ  S'n  ^ 


(B-7) 


{[Ro  *  Lc5."#  ]  -R„} 
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A  Nusselt  type  analysis  was  performed  during  laminar  film  condensation 
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between  the  condensing,  boiling,  sonic,  and  entrainment  limits  for  a  given 
heat  pipe  geometry. 
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